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i.  m        now 

The  maximum  power  level  of  a  nuclear  reactor  is  fre- 
quently limited  by  the  rate  of  heat  removal.   The  actual  de- 
sign of  a  power  reactor  represents  a  compromise  between  the 
requirements  of  heat  transfer  and  those  of  nuclear  physics. 
Nuclear  reactor  theory  is  not  sufficiently  accurate  to  t;uar- 
a,:tee  that  a  particular  reactor  (especially  one  of  a  new 
design)  will  £o  critical.   Therefore,  the  nuclear  characteris- 
tics of  a  reactor  design  must  be  obtained  before  construction 
of  the  full  scale  reactor.   Come  of  this  information  can  be 
obtained  from  the  use  of  a  subcritical  assembly. 

Experiments  with  uranium  graphite  subcritical  assemblies 
can  be  carried  out  to  determine  the  various  nuclear  constants 
of  a  proposed  reactor  lattic  .   These  lattice     ■tantl  are 
the  material  buckling,  multiplication  constant,  lattice  dif- 
fusion length,  and  thermal  utilization  of  the  unit  cell  in- 
cluding the  process  tube  assembly. 

In  this  investigation  six  lattice  configurations  were 
ildered,  and  the  buckling  of  each  was  determined  experi- 
mentally,  lattice  cells  of  6  in.,  0.5  in*,  and  12  i.u,  de- 
noted I,  II,  and  III  respectively,  constituted  the  thrcr  geo- 
metric lattice  arrangements.  For  each  of  these,  both  the 
?1wet!t  and  "dry"  cases  were  considered.   The  measurements  and 
calculations  with  water  present  in  the  cooling  annulus  were 


designated  as  "wet'1,  and  those  with  the  process  tube,  but 
without  water,  as  "dry". 

These  six  lattice  configurations  were  also  considered 
in  the  theoretical  analysis.   The  buckling  or  each  of  the 
six  lattice  configurations  was  determined  theoretically  by 
two  methods •   >nly  one  size  fuel  element  and  one  size  proc- 
ess  tubing  was  used  in  the  experimental  and  theoretical 
analyses  of  this  Investigation* 

The  methods  of  Murray  (1)  and  of  Humssy  and  Volkoff  (2), 
(3)  were  used  to  determine  theoretically  the  thermal  utiliza- 
tion and  the  resonance  escape  probability*   Murray1 a  method 
was  shown  to  give  less  conservative  results  tnan  the  method 
of  Fuasey  and  Volkoff,   The  theoretical  multiplication  con- 
stant derived  for  each  configuration  considered  was  based 
on  literature  and  calculated  values  of  the  various  diffusion 
lengths,  the  Fermi  age,  and  the  various  macroscopic  cross 
sections.  The  product  of  »7  and  £  was  assumed  to  be  es- 
sentially constant  Cor  the  various  configurations  investi- 
gated. 

Analyses  of  the  data  for  the  given  configurations  were 
carried  out  fcj  use  of  the  classical  formula  for  the  buckling 
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A  comparison  between  the  "wet"  end  "dry"  oases  of  the 
three  lattiee  cells  considered  was  presented  usin^  the  experi- 
mental an<*  theoretical  values  of  the  buckling*   By  an  analysis 
ol"  the  results  obtained  in  this  investigation,  the  effect  of 
the  water  coolant  on  the  nuclear  properties  of  the  various 
lattice  configurations  was  determined. 


II.   REVIEW  OF  LITERATI 

The  theoretical  calculation  of  the  multiplication  factor 
for  a  heterogeneous  reactor  with  coolant  is  an  extension  to 
the  method  of  calculating  the  thermal  utilization  and  reso- 
nance escape  probability  factors  for  an  infinite  heterogeneous 
uranium-graphite  reaotor  presented  by  Weinberg  (k)  »  15)* 
Guggenheim  and  Pryce  (6)  extended  the  method  to  include  a  non- 
absorbing  coolant  and  a  thin  sheath  of  absorbent  aluminum 
around  the  fuel  element.   Houston  (7)  extended  the  method  to 
Include  neutron  absorption  in  a  coolant.   Houston's  results 
are  more  suitable  for  £as  or  heavy  liquid  metal  cooled  re- 
actors due  to  his  assumptions  concerning  the  coolant's  negli- 
gible moderating  properties  and  small  volume  fraction  of  the 
lattice  cell.  Rumsey  and  Volkoff  (2),  (3)  extended  the  above 
methods  further  to  include  the  moderating  effect  of  a  water 
annulus . 

Murray  (1)  presented  a  simplified  method  of  calculating 
the  thermal  utilization  which  considers  all  but  the  uranium 
and  the  graphite  in  the  lattice  cell  as  poisons  which  reduce 
the  basic  uranium  graphite  thermal  utilization.    array's 
method  for  the  calculation  of  the  resonance  escape  probability 
does  not  consider  the  so  called  poisons.  The  Rumsey  and  Vol- 
koff method  appears  to  be  at  the  present  time  the  most  com- 
prehensive theoretical  analysis  of  a  water  cooled  uranium- 
graphite  lattice  cell. 
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The  experitaents  with  subcrltieal  assemblies  of  the  Han- 
ford  Atomic  Products  Operation  included  an  effort  to  determine 
the  effect  of  a  water  annulus  around  the  fuel  element  ra  the 
buckling  The  effect  of  aluminum,  air,  and  other  impurities 
in  the  lattice  upon  the  buckling  was  also  investigated  { 
A  few  measurements  were  taken  using:  internally  cooled  cylin- 
drical fuel  rods,  but  the  majority  of  the  coolant  simulation 
experiment!?!  were  with  solid  cylindrical  fuel  rods  of  natural 
uranium  surrounded  by  an  annulus  of  water*   The  objectives  of 
the  subcritical  assembly  experi'^nts  at  Kanford  (3)  were  to 
determine  exp+ylintally  the  various  lattice  constants.   The 
existing  lattice  theory  was  supplemented  and  improved  where 
possible,   lattice  Measurements  were  taken  with  three  differ- 
ent fuel  rod  diameters:  0.925  in«,  1.175  in,,  and  1,36  I  . 
The  lattice  spacing  and  the  water-aluminum-uraniua  ratios 
were  varied  Cov   each  slug  sisuu 

The  results  of  the  series  of  tests  using  both  "wet*  and 
"dry*  lattices  were  compiled  by  Clayton  (9).   A  comparison  of 
theory  and  experimental  results  for  one  fuel  element  diameter, 
Six  lattice  spacings,  and  the  "wet*  and  "dry"  confi^urationa 
for  the  Kanford  Atomic  Products  Operation  was  done  by  Gast 
and  others  (2).   Oast  used  the  method  of  Humsey  and  Volkoff  to 
calcalste  the  theoretical  values  of  the  thermal  utilization 
and  the  resonance  escape  probability. 

Related  conclusions  (8)  from  the     ord  work  with  water 
cooled  uranium  graphite  lattices  were: 


(a)  ?or  a  given  natural  uranium,  fuel  assembly, 
there  Is  a  Maximum  in  the  buckling  vs.  C/U 
atom  ratio  curve.  The  slue  diameter  and  the 
water  annulus  thickness  arc  inversely  pro- 
portional to  the  value  of  C/U  at  which  the 
maximum  occurs.  This  maximum  occurs  between 
C/TJ  values  of  50  and  100. 

(b)  The  buckling  "dry"  is  greater  than  the  buck- 
ling "wet"  for  larger  lattioe  spacing*  ( 
atom  ratio).  As  the  lattice  spacing  Is  de- 
creased, the  bucklin     ves  cross  and  the 
"wet"  buckling  becomes  greater  than  the  "dry" 
buckling.  For  this  lower  C/U  atom  ratio 
region,  externally  cooled  slugs  exhibit  a 
"fail-safe"  behavior  on  the  loss  of  water. 

The  following  conclusions  were  a  result  of  the  Hanford 
experl?nents(2) .  The   state  of  the  heterogeneous  graphite 
uranium  lattice  theory  at  present  is  not  completely  satis- 
factory. More  experimental  measurements  of  individual  par. 
eters  such  as  f  and  p  are  necessary  in  order  to  establish  the 
basis  for  an  improved  theoretical  analysis.   Corrections  for 
variations  in  reactor  material  densities  and  purities,  and 
allowances  for  holes,  such  as  for  control  rods,  need  further 
investigation  (2).  Finally,  the  treatment  of  the  energy  spec* 
trum  of  neutrons,  fast  to  thermal,  requires  more  accurate 
development  of  neutron  resonance  and  capture  reactions. 
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A.   General  Procedure 


fhe  product  of  the  fast  fission  factor,  €  ,  the  primary 
fission  neutron  factor,  t)    ,  the  thermal  utilization  factor, 
t ,  and  the  resonance  escape  factor,  p,  ia  the  infinite  multi- 
plication constant,  kv   . 


6V|fp  *  k„ 


(2) 


A  theoretical  analyaia  of  eac     the  factors  would  have 
been  extremely  complicated  without  certain  simplifying  as- 
sumptions.  The&e  assumptions  are  stated  where  applicable  in 
the  following  theory.   Each  of  the  factora  was  considered 
separately,   fhe  product  €V|  was  determined  and  assumed  con- 
stant for  all  the  lattice  configurations  considered,   repres- 
sions for  f  and  p  were  determined  using  the   methods  of  Hurray, 
denoted  Method  A,  and  the  methods  of  Bumsey  and  Volkoff ,  de- 
noted Method  Bf  for  the  "wetn  and  "dry*  cases  of  Lattices  I, 
II  and  III. 

Using  values  from  the  literature  for  the  lattice  diffu- 
sion length  and  the  Fermi  age  with  the  theoretically  deter- 
mined infinite  multiplication  constant,  the  buckling  was  deter- 
mined fro™  KQuation  1.   The  theoretical  buckling  was  determined 
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for  lattices  I,  II,  and  III  for  both  the  "wet"  and  "dry" 
cases  using  both  Methods  A  and  B, 

The  effect  of  the  air  annulus  and  air  holes  on  the  neu- 
tron density  in  all  lattice  jeome tries  was  neglected  due  to 
the  magnitude  of  the  effect  in  comparison  with  other  material 
effects  in  the  lattice. 

Subscripts  used  in  this  analysis  were  defined  as  follows: 

u   natural  uranium 

al  slug  can,  25  aluminum 

w   water 

g   graphite 

m   moderator 

p   process  tube,  61S  aluminum 

a   absorption 

s   scattering 

re  resonance 

sa  scattering  "absorption" 

f   slow  neutron  fission 

t   transport  theory 
Any  deviation  of  notation  from  the  above  list  was  individual- 
ly defined, 

B,  Past  Fission  'actor 

The  fast  fission  factor,  €  ,  is  the  ratio  of  the  total 
number  of  neutrons  produced  by  fission  to  the  number  of  neutrons 


' 
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produced  by  fission  by  thermal  neutrons.   Assuming  that  the 
primary  neutron  souree  is  distributed  approximately  uniformly 
over  the  fuel,  a  simplified  version  of  the  equation  ^iven  by 
Glasstoae  and  Sdlund  (5,  ?•  278)  for  the  fast  fission  factor 
is 


<a. 


e  •  x  •  -*-  ,  ..^:  \  l.  *  ■ ,   P   , ,    *  _,  (3) 


where  V       is  the  probability  that  a  fission  neutron  born  in 

a  rod  will  make  a  collision  inside  the  rod  in  which 
it  was  created! 
<Ta  Is  the  absorption  cross  section  for  fast  neutrons, 

that  is,  the  sum  of  the  fission  and  the  capture  cross 
sections j 
g-   is  the  elastic  cross  section  for  fast  neutrons; 
(J\   is  the  inelastic  cross  section  for  fast  neutrons j 
and   (j"t  is  the  total  cross  section  for  fast  neutrons. 

<rt  -  <ra  ♦  <re  *  cr.. 

Using  fast  neutron  cross  section  values  as  given  by  Gug- 
genheim and  Tryce  (6,  p.  5D,  Equation  3  reduces  to 


e  •  x  «  °-°W3  '■'  .  (h) 

1        1  -  0.521  W 


! 


■ 


V 


to 

Those      leal  coefficients  are  such  that  when  the  fuel  ele- 
ment radius  is  less  than  1»7  cm  the  calculated  value  of  €- 
agrees  reasonably  well  with  experimental  results  obtained  by 

onne  National  Laboratory  (6,  p*£l)»   V  or  natural  uranium 
of  density  19  g/ci?r  and  a  rod  radius  oC   &*$  in,,  the  value  of 

*<n.  r.  is  0*262,  where  ?f  is  the  number  of  uranium  atom*  per 
%     xx, 

unit  volume,  and  r   ia  the  radius  of  the  fuel  rod.   The  quan- 
tity P  is  a  iunction  of  §(T#  r  »  From  a  relationship  of  i  and 

v  u 

HO",  r,  shown  by  Guggenheim  and  Pryce  (6,  p*  $1)   and  the  above 

v   U 

conputed  value  of  HO"t  r  ,  the  value  of  P  determined  was  0.25* 

Solving  Fnufition  k   for  €  gave  a  value  of  1.027*   This 
value  for  the  fast  fission  factor  was  assumed  constant  for  all 
configurations  investigated* 

C*   Primary  Fission  Neutrons 

The  number  of  fission  neutrons  released  per  thermal  neu- 
tron captured  in  the  fuel  is  the  factor,  n  •   Glasstone  and 
Edlund  define  h   (5,  p.  B3)   as 

u 


where  V     Is  the  average  number  of  neutrons  produced  per  thermal 
neutron  fission; 
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and  £   is  the  total  cross  section  for  absorption  of  thermal 


•  u 
•a. 

neutrons • 


The  microscopic  thermal  fission  cross  section  for  natural 
uranium  is  i^lS  barns.   Applying  the  raaxwelllan  distribution 
correction  of  1.128  and  the  "not  -  l/v"  factor  of   »ft!  for 
XT-*-*   to  this  value  gives  an  effective  fission  cross  section  of 
3.635  barns  for  natural  uranium.   The  microscopic  thermal  ab- 
sorption cross  section  for  natural  uranium  is  7.68  barns.   Ap- 
plying the  maxwell ian  distribution  correction  and  a  "not  -  l/v" 
factor  of  0.990  for  natural  uranium  to  this  value  determines 
an  effective  absorption  cross  section  of  6.7k  barns  for  natu- 
ral uranium  (1,  p.  32). 

Murray  (1,  p.  98)  gives  a  value  of  2.1+6  for  P    .   Sub- 
stituting the  above  quantities  in  Equation  $,   the  calculated 
value  of  r)    was  1.327. 

It  is  normal  practice  in  working  with  subcritical  assem- 
blies to  experimentally  determine  kw  and  then  divide  k^  by 
the  product  £  pf  to  arrive  at  a  value  of  in  .   The  value  of  in 
found  by  subcritical  assembly  experiments  was  1.3 0d  neutrons 
per  thermal  neutron  captured  (10,  p.  85) •   This  value  is  in 
reasonable  agreement  to  the  computed  value  of  1.327.   The  val- 
ue of  iri  used  in  the  following  analysis  was  1.3)8  and  was  as- 
sumed constant  for  all  configurations  investigated. 


r« 


12 
D.   Thermal  utilisation 

In  heterogeneous  lattice  theory  the  finite  lattice  re- 
actor is  assumed  to  be  an  infinite  lattice  array.   The  prob- 
lem of  calculating  t   and  p  is  simplified  b  replacing,  the 
square  cells  by  equivalent  cylindrical  unit  cells.   Table  1 
gives  the  equivalent  cell  dimensions  for  Lattices  I,  II,  and 
III.   Assure    -hat  little  tfe*Qf«  in  the  overall  neutron 
flux  is  experienced  in  traversing  one  cell,  an  individual 
cell  nay  be  representative  of  the  total. 

The  methods  I    i^ray  (1)  and  or  }  umsey  and  Volkoff  (2), 
(3)  were  used  to  determine  theoretically  the  thermal  utiliza- 
tion.  These  methods  were  considered  separately  and  were  de- 
noted Ketho4a  k   and  E  respective!  . 

1. 

in 

The  depression  of  the  thermal  neutron  flux  in  the  fuel 
and  adjacent  moderator  complicates  the  calculation  of  the 
thermal  utilisation.  Average  flux  and  volume  weighting  fac- 
tors must  both  be  included  in  the  fractional  absorption  of 
the  fuel.   From  the  spatial  variation  of  the  flux  In  a  cell, 
the  average  values  of  the  flux  can  be  determined. 

The  ther-nal  neutron  diffusion  equations  for  the  uranium 
and  the  moderator  from  diffusion  theory  are 


K*  K-K  zu  +  su  =  °  <6> 


13 
and 

D«  V2    K    *    K      ^m    *     St,  "  °  (7) 

m  v       'm    w   m    ni 

Equations  6  and  7  can  b©  simplified  with  th©  following  as- 
sumptions: 

(a)  The  variation  of  neutron  flux  at  a  given 
eel!  radius  will  be  the  same  everywhere. 

(b)  There  are  no  thermal  neutrons  produced 
in  the  uranium,   3  is  zero. 

(c)  The  production  rat©  of  thermal  neutrons 
in  the  moderator  is  independent  of 

position.   3  is  constant. 

m 

(d)  The  flux  is  assumed  not  to  vary  along 
the  fuel  cell  axis. 

Using  the  above  assumptions  and  cylindrical  coordinates, 
Equation  6  for  the  fuel  becomes 


^2    x  \  dx  / 


-  A,  -  0  (8) 


where  x  ■  $  ^r  and  tf   is  the  inverse  diffusion  length  of 
uranium,  and  Equation  7  for  the  moderator  becomes 


x^MdrV    -  't.X—r?  «9) 


Ik 


where     D     la   the   diffusion  coefficient   for    the  moderator. 

rs 

m 

*f   is  the  square  of  the  inverse  diffusion  length, 
in 

Equation  0  for  the  fuel  is  the  modified  Bessel  equation 

of  zero  order.  T qua t ion  9  for  the  moderator  is  the  inhomo- 

geneous  type  of  the  modified  Bessel  equation  of  zero  order# 
The  solutions  of  Equations  8  and  9  are 

j*u(r)  =AI0  (tfur)  (10) 

and      ftjr)   *  C  Kn  («  p)  +  C  L  (  tf  r)  •  *■        (11) 

*-m 

respectively,   The  boundary  conditions  used  in  the  evaluation 
of  the  constants  a,   ,  MMl  fl  are  neutron  flux  and  current 
continuity  at  r  »  r   (uranium  rod  radius)  and  the  condition 
that  dtf/dr  «  0  at  r  *  r2  (the  equivalent  outer  radius  of  the 
graphite  coll). 

Applying  the  condition  of  zero  neutron  current  at  r2  in 
Equation  11 


h  <*m*y 


0  «  c  t^-7 — £*.  (12) 


thus  the  moderator  flux,   r.quation  11,   becomes 

1 
j*m(r)   m  g  m0   (jf^)    ♦  Ji  (13) 

m 


% 
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where 


«0  <**r)  *  *0  <*mr)  f  i;  UVL)      l,   (^mr) 

l    m  d. 


Continuity  of  flux  and  current  st  ru     |  &  pair  of  equa- 
tion* from  Equations  10  and  13  which  can  be  solved  simul- 
taneously to  give 


ft: 


m  **  m  1    ra  u 


(U) 


and 


s  — 


\   *u  Xl  (*uru> 


(15) 


where  K,  (^  pJ 

1    ra  2 


and 


m 


■ 


"•  *0<  "uru>  V  rf«ru>  +  r'u*u  VW  K)(^Pu» 


with 


dKn<x) 


The  thermal  utilization,  ff  is  the  ratio  of  the  thermal 
neutrons  absorbed  in  the  fuel  to  the  total  thermal  neutrons 
absorbed.  This  ratio  is  equivalent  to  the  ratio  of  the  ther- 
mal neutrons  absorbed  la  the  fuel  to  the  total  thermal  neutron 
source  assuming  no  net  leakage  from  the  cell# 


%i 


le 


u 


0 


/u  ^u  2»  r  <*r 


V   3 
m  in 


(16) 


fining  I  qua t ion  10  and  16  Into 


r  =» 


/  JL.\  2n u 

I  3  /    V 


u 


/   IQ  («^ur)  r  dr 


(17) 


and  3olvlng  the  Integral  (1,  p.  305)  £ives 


2*  Z 


f  ■ 


A  \  2*  2u    yu  *i  <  W 


(16) 


Combining  Equations  II4.  and  18  and  solving  Tor  the  reciprocal 
of  f  gives 


1 
? 


»  1  + 


v   r 

m ra_ 

u  *-u 


+  (E  -  1) 


(19) 


where 


F  ■ 


Vu 


I„  (  #   r  ) 

Q    ^  u. 

Xl  <  W 


(20) 


aid 


E  -  1  « 


Til 

V" 

u 


^  r 

m  u 


-  1 


(21) 
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,  is  the  disadvantage  factor  of  the  uranium  which  ia  th© 
ratio  or  t  e  neutron  flux  at  the  rod  surface  to  the  average 
flux  in  the  rod  interior  (5,  p«27   . 


The  average  thermal  neutron  flux  in  the  moderator  is  greater 
than  the  thermal  neutron  flux  at  the  rod  surface  because  the 
diffusion  coefficient  in  the  moderator  is  finite*   The  ad- 
ditional absorption  is  measured  b?  the  quantity  {E  ~  1),  the 
excess  absorption  term. 

Hurray  (1,  p*  100)  assumed  that  any  poisons  that  are 
tolerable  do  not  appreciably  disturb  the  basle  neutron  flux 
distribution  from  that  of  the  com  1  oration  with  fuel  and 
moderator  only.   Therefore,  the  aluminum  cladding,  coolant, 
and  the  process  tube  assembly  can  all  be  treated  as  poisons* 
The  reciprocal  thermal  utilization  with  this  modification  can 
be  written 

s_  poisons 

L_     u  **u  ^u 

1=0 

where  the  average  neutron  flux  in  the  added  poison  or  absorber 
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Assuming  that  for  the  cladding,  water  annulus,  and  proc- 
ess tube  the  ratio  J^/  #u  is  approximately  equal  to 
P     ir  )/  0  ,   Equation  25   can  be  solved  directly.   Kurray  (1, 
p.  100)  suggested  that  a  slight  refinement  to  1 ouation  23  may 
be  made  by  letting  Z.  be  the  difference  between  the  actual 
poison  cross  section  and  the  moderator  cross  section  that  it 
physically  replaces. 

Littler  and  Raffle  (10,  p.  9U)  §tn  a  useful  approxima- 
tion for  the  excess  absorption  term,  E  -  1. 


-  1  s  1'.   ^ 


r2  ~  rl 


In 


—is  .  i  -f 


2 


Ijr, 


(2W 


The  constants  used  in  T.  qua  t  ion  23  were  listed  in  Tables  1 
and  2.   The  equivalent  cell  dimensions  were  derived  from  di- 
rect measurements  and  calculations.   The  thermal  neutron  macro- 
scopic cross  sections  for  aluminum  (cladding),  water,  and 
graphite  were  calculated  from  Equation  25  using  literature 

values  of  CT   (1),  (5),  and  (11). 

n 


5   g  °* ! 


(25) 


where  Wm   is  6.02  x  1023  nuclei/mole, 


and 


A  is  the  mass  number. 
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Zj  for  the  aluminum  proceas  tube  was  calculated  ItlN  a 
modification  or  £quation  25  to  account  for  the  alloying  ma- 
terials in  the  aluminum  alloy.  The  raaxwellian  distribution 
correction  and  a  "not-l/vrt  correction  was  used  to  determine 
an  effective  absorption  cross  section  i'or  natural  uranium. 
All  the  macroscopic  cross  section  values  used  in  this  theo- 
retical analysis  are  shown  in  Table  2. 

Due  to  the  limitations  of  simple  fllffusion  tneory,  trans* 
port  theory  raust  be  used  for  the  calculation  of  Jf  .      Murray 
(lt  p.  88)  gives  for  the  uranium  fuel  the  relationship 


■J  *  tanh  —§  (26) 


whore  £  *8  fcke  total  cross  section, 

Solving  Equation  2$   usias  cross  section  values  from  Ta- 
ble  2   grave  a  value  of  0.675  cm   for  if   .   The  inverse  dif- 
fusion lengths  for  graphite  and  water  were  found  usin,;  I  qua- 
tion  27  and  values  of  the  diffusion  lengths  from  the  literature 
(5). 


*•$«  {ZUMX/Z 


un 


The  thertaal  utilizations  for  the  dry  lattice  configura- 
tions were  calculated  by  omitting  the  term  for  water  in  Equa- 
tion 23*     The  volume  ratios  in  Ecuation  19  and  2}   are  equiva- 
lent to  area  ratios. 


ri 


•  ^ 
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The  constants  listed  In  Tables  1  and  2  were  need  in 
Equations  19,  2%  23,  and  Zk   to  evaluate  the  composite  ther- 
mal utilizations  fop  the  various  lattice  co     rations  con- 
sidered.  A  Bessel  function  table  (12)  was  used  to  evaluate 
the  disadvantage  factor,  f9   in   I  qua t ion  20*  Equation  2k   was 
used  to  evaluate  the  excess  absorption  term,  (7  -  1). 

The  resulting  thermal  utilizations  computed  by  Method  A 
are  listed  in  Table  3  Tor  the  various  lattice  configurations. 
2.   Method  S 

A  second  method  of  calculating  the  thermal  utilization 
of  uranium  is  the  met  od  of  Rumsey  and  Volkoff  (2),  (3). 
method  Is  an  extension  of  prior  work  in  heterogeneous  lattice 
theory  (k) ,  (6)  and  includes  the  moderating  effect  of  a  water 
annulus  in  the  lattice  cell. 


fu«f   (1*6) 


(28) 


where     4--l»K.*R  +  R   ♦  R   ♦  S  +3  «■  B   + 

TZ  al    p    w    g    w    g    wp    we 


(29) 
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•  "relative  absorption"  terra,  R.,  denotes  per  thermal 
neutron  captured  In  the  uranium,  the  number  of  thermal  neu- 
trons which  would  be  captured  In  the  jth  medium  if  the  neutron 
density  In  the  jth  medium  were  uniformly  equal  to  the  neutron 
density  at  the  uranium-aluminum  interface. 


u  w  u  u 


where  F  is  the  disadvantage  factor  of  the  uranium,  Equation 
20* 

The  "self  blocking"  term,  S-,  denotes  per  thermal  neutron 
absorbed  in  the  uranium,  the  excess  number  of  thermal  neutrons 
captured  in  the  Jth  medium  due  to  the  exeeas  density  of  neu- 
trons in  the  jth  medium  over  the  neutron  density  &t  the  i  -  jth 
interface* 


3w   2   w  K  OZ) 

and      3g  -  (1  -  1)  [l  ♦  Bftl  #  Hp  ♦  Rw  ♦  B^  ♦  3w]   (33) 

where  (E  -  1)  is  the  disadvantage  factor  of  the  moderator  or 
excess  absorption  term  given  by  Equation  2h* 

The  "blocking  term",  B..,  denotes  per   thermal  neutron  ab- 
sorbed  in  uranium,  the  excess  number  of  thermal  neutrons 
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captured  in  the   Jth  medium  due  to  neutron  density  rise  across 
the  ith  medium* 


w 


<3W 


B         «      «*f    2    t   ; 
W£  ^  W         W 
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(35) 


In  the  abov*  six  eouations  #     ia  the  inverse  diffusion 
length  of  water,  and  t  is  the  thickness  of  the  water  annulus. 
The  cross  section  valuer  ara  listed  in  Table  8t  and  the  volume 
values  per  slug  are  listed  in  Table  1.  The  value  of  (  -1) 
is  given  by  Equation  2k  which  is  an  approximation  of  the  M  - 
plex  Bessel  function  in  quation  21,   The  approximation  is 
considered  to  be  sufficiently  accurate  for  the  two  lattj 
geometries. 

rhe  production  rate  of  thermal  neutrons  per  unit  volume 
per   second  for  water  and  graphite  are  denoted  o  an     ie- 
apectively.  The  ratio  of  the  slowing  down  power  of  water  to 
graphite,  Q^At  »  equals  20  {3,  p«  22), 

igure  1  t»d     are  2  show     leutron  density  distribu- 
tion in  the  dry  and  wet  lattice  configurations  respectively. 
The  several  terms  of  the  cotnpetitlve  absorption,  Equation  29, 
are  shown  schematically  in  these  diagrams.  The  placement 


rtqao 


the  individual  absorption  terms,  defined  in  Equations  31 
through  35*  and  the  relative  magnitude  of  these  various  terms 
©an  also  be  seen  in  the  diagram.  Any  absorption  of  neutrons 
by  the  air  in  the  air  annulus  was  neglected.   The  net  flow  of 
thermal  neutrons  across  the  boundary  of  a  cell  was  assumed  to 
be  »ero« 

The  terras  B^,  Sp,  B^,,,  B.lp,  i^,  and  Bpg  which  are 
shown  in     ire  1  and  Figure  2  were  assumed  a::       ible  and 
were  not  included  in  Equation  29 •  At  the  uranium-aluminum 
interface,  continuity  of  the  neutron  flux  density  and  current 
density  was  assumed*   The  neutron  current  density  was  also  as* 
sumed  as  linear  through  the  thin  mediums  (al,  ;;,  or  w)  con- 
sidered. 

As  an  example,  in  deriving  the  I  ,  term,  the  aluminum 
surface  neutron  flux  is 


I.™-  .  ^(ru,  .  t.x 


\-3r7 


r 

u 


and  with  Equation    (10)    is 

By  definition  the   *•#!£  hloekinr:"  tern     for  aluminum  ia 


/ 


( hl  &ve        m  eurface)  ■  r-  al 

3      „      l'*l       'A L      V»l    \  (36) 

AT6         *u  c 


and  taking  from  (10,  p.  89) 


and  substituting  gives 


S.l£|*al2  *«12  «?> 


Similarly 

1  £  . 

w   2  w    w 


|  £  1  t  2  *  2  (32) 


and 


3  Sit2  <*  2 

?   *  P    P 


2  . 


A  calculation  of  S  - ,  I  qua t ion  37,  showed  that  with  tf   , 

0.OO2UI4.  en!*2  (X'0f  p.  93)  and  t  x  *  0.102  en,  that 

S  -  «  1.27  x  10  .  This  value  is  of  a  magriituds  which  can  be 

considered  negligible* 

a 
The  very  small  factor,  (  Jf.  t,)1",  Tor  the  cladding  and 

process  tube  causes  the  ter?ns  9L%$        .  ^aiw#  ^l<  »  B  .  ,,  and 

B   to  become  negligible  with  respect  to  the  other  terms  in 
P£ 

■ure  1  and  .-'igure  2.   With  these  eliminations,     -e  1  and 
Figure  2  agree  with  Equation  2  . 
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The  constants  lis  to*4!  In  Tables  1  and  2  were  used  in 
Equations  2&   through  3£  to  avalua     e  wet  configuration 
thermal  utilizations  by  the  Rumaey  and  Volkoff  method.   cus- 
tion  20  was  used  to  evaluate  ?  in  liquation  31.   The  value 
(   -  1)  in  Fquation  3^   was  evaluated  by  using  Fquation  2I+. 
The  thecal  utilizations  for  the  dry  configurations  were  cal- 
culated using  Equations  28  through  35  with  the  terms  for  wa- 
ter omit   r  ' .      The  resulting  I    ftl  utilizations,  Method 
are  listed  in  Table  3  -'or  the  various  configurations, 

E.   Resonance  Escape  Probability 

The  fraction  of  fast  fission  neutrons  that  reach  thermal 
energy  without  experiencing  capture  while  slow     own,  is 
oalled  the  resonance  escape  probability,  p.   The  methods 
Murray  (1)  and  of  r umsey  and  Volkoff  (2)  were  also  used  to 
theoretically  determine  the  resonance  escape  factor .   These 
methods  are  considered  separately  and  are  denoted  Methods  A 
and  B  respectively. 
1.   Method  A 

Murray  (1,  p.  92)  ^ives  the  hetero-  eneous  equivalent  for 
p  as 

P  ~  **P       -  a  U  ./**  - 
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where  21^  is  the  effective  macroscopic  cross  section. 

The  form  of  Venation  38  suggests  a  resemblance  to  similar 
ratios  appearing  In  the  expression  to  evaluate  f. 

A  definition  of  the  resonance  utilization,  t    ,  Is 

*     a resonance*  absorption  In  fuel 


r   resonance  absorption  in  fuel  ♦  scattering  "absorption" 

by  moderator 

(J 

The  resonance  absorption  by  the  uranium  Is  ft     V  "Z. 
from  a  the-mal  flux  absorption  a  at.  I  \,  ,  tihe:   £  is  an  ap- 
propriately chosen  average  cross  section  of  uranium  over  the 
entire  resonance  flux  region. 

Assuming  that  the  flux  over  the  resonance  r    a   Is  l/B 
dependent,  the  average  croas  section  is 


u 


U  dE 
J 


/cra'J  -    o-p 


In  ~ 


1 
where  In  r—  «  5»6  for  uranium  metal  (5#  p»  2?3).   The  presr 

"2 
of  the  moderator  and  the  geometric  natu      the  heterogeneous 

lattice  requires  that  an  effective  resonance  integral  (0~  ) 

re 

be  used  for  the  value  o£'G"r-      I  qua t ion  1+0  is   then 

re 

0"u  =  TX  <**> 


"7 
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where  (1,  ?•  93) 

(j  ■  m  /   (  (T  )  eff  ~  «  9.25  •   i   /«       di2) 

re        w 

and  s/m  Is  the  fuel  surface  to  mass  ratio  in  on  /  . 

o~  u    zT  u 

fhe  scattering  "absorption"  by  the  moderator  would  be 
Hi  *■   Z»     from  a  therml  flux  absorption  analo  ,  , 


where  <r-a      y    *   *-  8  K  „.      *   *- s  /,.•  » 


With  quations  39,  U3,  and  I4I4.  the  magnitude  of  the  exponent 
of  p  In  Equation  38  is 


*»  **a 


M" 


u  *u  *re  (J.5} 


Murray  (1,  p.  9U)  gives  an  empirical  expression  for   the  ef- 
festive  resonance  inverse  diffusion  length  of  uranium, 


tf£C  -  0.0222^     cm"1  (1^) 


where  p     is  expressed  In  g/cisr  . 
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Th©  effective  resonance  inverse  diffusion  length  for 
graphite  is, 


re 


I 


«T"i 


*   ,  V;iZ  7) 


It 


where  Z.  is  the  actual  transport  cross  section  of      r&ph- 
t 

its  averaged  over  the  entire  resonance  flux  re 
evaluate  £..,  the  relationship 


Z<.  m    I"  (1  -  cos  Q  )  V 

Is  s 


is  derived  from  (5,  p»  98) 

*4-  s  KA  m   cos  %  (k$) 

where     cos  '#  ■  2/3 A  0) 

and  A  is  the  moderator  mass  nusaber. 

The  method  of  calculating  f  is  now  borrowed  to  find  f  • 

r 

ation  19  modified  is 

i-  *  1  *  vfi-  ■  ■  ■-  P  ♦  (   -  1 )  1) 

xr       u  *0 

The  values  of  *™   and  tfr*  were  determined  from  Equations  1* 

u        | 

and  lf7.   2",,  Mi  Z2!1  were  calculated  from  Equations  k3   and  Mt* 
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Using  the  above  values  of  Jf r     and  #*   ,  P  was  determined 
from  Equation  20  and  (i  -  1)  was  determined  from  Fouation  2l|. 
Table  2  lists  all  the  important  lattice  material  constants. 
Solving  Knuation  £l  and  using  the  value  of  f   in  T qua- 

8 

tion  hS   determined  the  exponent  of  p  in  Equation  33»  The  res- 
onance escape  probabilities  for  the  various  configurations, 
Method  A,  are  listed  in  Table  3*   The  values  of  p  for  the  wet 
and  dry  configurations  were  the  same  using  Method  A. 
2.   Method  B 

Volkoff  and  Kurasey  (2,  p.  293)  express  the  equation  for 
the  resonance  escape  probability  as 


exp 


fr 


1  •  f ' 

r 


u 


(52) 


where  f  is  the  resonance  utilization  of  uranium  for  reso- 
r 

nance  neutrons.   Noting  the  form  similarity  of  the  exponent 
of  Equation  52  to  the  thermal  utilization  equations,  the  fac- 
tor, f  ,  is 
r 

fr  «  fQ  (1  >  5)  (53) 


with 
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ohanging  the  cross  sections  7  to  2^  and  Z.^  to  Z-  In 

nations  31,  33,  and  35   gives  a  solution  to      ion  5k> 
6  is  riven  by  Equation  30. 

Substituting  the  individual  terms  in  Equation  5k   gi*M 
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where 
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Combining  Equations  5.2,  53,  and  55  rives 


p  35  exp 


.   (1  +  6) 

Z  +  Y  - 


(53) 


The  value  of  I  was  again  obtained  from  Equation  2ij..   The 
average  cross  section  of  uranium  over  the  resonance  region, 
Z0,  is  given  by  Equation  U3 •  The  scattering  "absorption" 
cross  section  of  the  moderator  over  the  resonance  region,  Z.-., 
is  given  by  Equation  IJ  .   Using:  values  from  Tables  1  and  2 
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and  Equations  I     »  >©,  "3,  ;,  ,  l-/->   I  ft    $4,  I  7»  and  58  the 

resonance  escape  probability  was  determined  for  Method  I  . 
These  values  of  p  are  listed  in  Table  3  for  the  lattice  con- 
figurations considered.   The  dry  values  for  p  were  obtained 
by  eliminating  the  terms  for  water  in  all  equations, 

.   Buck 11 

The   theoretical  buckling   was  computed  using  Equations  1 
and  2,      The   diffuaio  i   Tor   each  lattice  BO  ration 

was  calculated  by    (5,   p» 


t.2  «  %-  -  L  2(1  -  f )  (59) 


where  r.  is  the  diffusion  length  in  pure  moderator ♦   The  val- 

I 

ue  of  f  used  in  each  case  was  the  theoretically  determined 
value  listed  in  Table  3.   The  Fermi  age,  r  ,   of  the  thermal 
neutrons  for  the  lattice  configurations  was  assumed  to  be  the 

same  as  t   for  the  moderator.   Murray  (1,  p.  123)  gives  the 

p 
value  of  *>*  as  361|  cm  .   The  value  or  kw  for  each  e     na- 
tion was  determined  from  the  values  o-  €  ,  Y)    ,  f ,  and  p  listed 
in  Table  3  for  each  con: i  .oration. 

The  volume  of  a  cubical  critical  reactor  having  the  | iven 
lattice  co     nation  is 

V-  ■  2^\  A 

g 
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The  square  diffusion  length  and  the  buc 

W   e&c     ;.he  lattice  configurations  considered  are  listed 
in  Table  !}•    o  volume  and  cube  side  length  of  a      al 
critical  reac  const      '  ittice  configurati 

are  also  listed  in  Table  iu 
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iv.    nmn  igatioi 

A.   Description  of  Equipment 

A  picture  of  the  subcritical  assembly  with  a  8*$   in. 
lattice  spacing  is»  shown  in   i.ure  3.    e  physical  dimen- 
sions of  the  graphite  assembly  were  60  in.  x  60  in.  x  79  in. 
The  actual  graphite  used  in  the  construction  of  the  assembly 
was  originally  7  in.  diameter  solid  praphite  rods.   These 
rods  were  squared  off  to  a  6  in.  x  6  in.  cross  section  (3.5 
in.  radius  on  the  corners)  and  stacked  as  shown  in  Figure  3« 
The  top  five  rows  of  the  graphite  rods  in  the  assembly  were 
5  in.  x  6  in.,  and  tne  bottom  nine  rows  were  6  in.  x  6  in. 
■  density  of  the  graphite  was  determined  from  a  wei 

and  dimension  measurement  of  a  mass  of  graphite.   The  graph- 

3 
ite  density  was  found  to  be  1.5>6  f/iri 

The  wraphlte  assembly  was  covered  on  the  top  and  four 
sides  with  a  10  mil  sheet  of  cadmium  sandwiched  between  a 
.375  in.  thickness  of  plywood  and   0.125  in*  of  masonite. 
The  cadmium  sheet  gave  an  effective  "black  wal  si 

neutrons.   &4        .2   in.  total  spacing  allowance  for  the 
graphite  rods  a  horizontal  section  between  the  "black  walls" 
of  the  assembly  was  60«5  i««  x  t %$   lru  in  area.   The  assem- 
bly was  mounted  on  a  concrete  base,  and  tanks  of  water  were 
placed  under  the  graphite  structure  to  insure  that  adequate 


shielding  was  effected  with  the  neutron  30urce  In  place. 
Figure  3  shows  the  relative  placement  or  these  tanks* 

The  lattice  spacing  was  varied  by  changing  the  •     ©- 
ment  of  the  uranium  filled  process  tubes  in  the  assembly. 
The  6  in.  lattice,  Lattice  I,  waa  the  oo.      ation  with  all 
the  lattice  holes  filled  with  process  tubes.   The  6«$  In. 
lattice,  Lattice  II,  waa  the  configuration  as  shown  on 
ure  3.  By  re  ■■■novin<j.  alternate  rows  of  the  process  tubes  from 
the  8.5  in.  lattice,  the  12  in.  lattice,  Lattice  III,  was 
realised. 

lve  plutonium-berylllum  neutron  sources  contained  in 
individual  ridit  cylindrical  containers  were  placed  under  the 
center  of  the  graphite  assembly.   Each  of  these-  sources  emit- 

z. 

ted  approximately  1.63  x  10  neutrons  per  second;  therefore 
the  total  source  strength  was  approximately  8.l£  x  10  neu- 
trons per  second.  A  cruciform  fo  one  try  was  chosen  for  the 
placement  of  the  five  sources.  One  source  was  at  each  of  the 
cruciform  ends,  and  one  waa  placed  at  the  cruciform  center. 
The  total  source  was  in  this  way  effectively  contained  in  a 
circle  of  about  3.5  in.  in  diameter. 

The  cylindrical  uranium  slugs  were  I.OSO  in.  in  diameter 
and  8.l}.0  In.  In  length.   Assuming  a     11  cladding  and  a 
200  mil  slug  end  cap  of  23  aluminum,  the  actual  uranium  fuel 
size  was  a  rod  1.    .  .  .   ianeter  and  ..  0  in.  in  length. 
density  of  the  uranium  was  de  rom  the  weights 
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and  dimensions  or  several  random  sample  el    .   ?he  density 
•f  natural  uranium  determined  in  this  way  was  19.   /«r« 

The  smallest  lattloe  spacing  investigated,  the  6  i  . 
lattice,  required  117  sections  of  6061-76  (commercial  desi.  - 
nation  61S)  aluminum  tubin  .   The  tub.       t-ions  had  an 
outer  diameter  of  1.375  in,,  a  wall  thickness  of  35  mils, 
and  a  length  of  62  in.   The  compos! tic      ]      liana  is 
0*2$%   cop;  er,  0.60^  silico •  ,  l.QO£  magnesium,  V.25%   chromic 
urn,  and  the  remainder  aluminum. 

The  spacing  wire  composition  was  23  aluminum  and  had  a 
0*102  in.  diameter.   The  wire  was  used  to  center  the  uranium 

s  in  the  process  tube*   Seven  uranium  slugs  were  loaded 
in  each  tube  with  a  helical  wrapping  of  wire  around  the  slugs 
for  spacing  in  the  tube*   Approximately  10  feet  of       ,'as 
used  for  each  seven  slu,  . 

Figure  if  show©  a  scale  drawing  of  an  actual  unit  cell 
of  the  subcritical  assembly  with  a  loaded  process  tube  in 
place.   -lumber  seven  rubber*  stoppers  were  used  at  both  ends 
of  the  process  tubing  to  contain  the  water  coolant.   The 
stoppers  were  also  used  for  the  dry  configurations  to  help 
contain  the  uranium  slugs  during  ha;      of  the  loaded  proc- 
ess tubes.   Allowing  two  process  tube  lengths  of  wire  for 
each  tube,  the  effective  thickness  of  the  process  tube  was 
calculated  to  be  1;.0  mils.         active  thickness  of  the 
water  annulus  was  then  found  to  be  ")*107  in. 


36 

Indium  foils  were  mounted  on  :'lat  aluminum  plates  which 
in  turn  were  placed  on  an  aluminum  holde*  .     •  holder  was 
designed  so  that  the  depth  of  penetration  of  each  Toil  into 

e  slots  in  the  assembly  could  be  accurately  determined. 
The  indium  foils  were  1.5  in.  x  1.0  in.  x   .   3  In. 
en  average  weight  of  C.59^3  grams.   The  average  thickness 
of  these  foils  was  6.1lv  mg/em  .  Figure  5  shows  the  foil 
holder  and  foil  placement  during  count    .       eometry  of 
the  foil  with  respect  to  the  rlass  wall  counting  tube  was 
held  constant  by  use  of  the  counter  shelf  holder  shown  In 
-ure  5. 

A  model  181A  scaler  manufactured  by  Suclear-C     •  and 
e  glass  wall  Gleger  tube  were  used  to  count  the  beta  activity 
of  the  irradiated  indium  foils.   A  conventional  stop  watch 
was  used  for  the  tir     r  the  counting  period. 

B.   Procedure 

The  uranium  filled  process  tubes  were  loaded  into  the 
graphite  assembly  In  varying  arrangements  for  the  three  geo- 
metric lattices  invest!  ated.     I  annulus  between  the  slug 
and  the  tube  wall  was  filled  with  water  for  the  wet  c      ira- 
tions. 

The  neutron  tlux  distribution  was  determined  by  the  indi- 
um foil  activation  method.   L-y  determining  the  activity  of  a 
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foil  after  a  lonr  p«     of  Irradiation  tad       an  al- 
lowance for  decay  during  the  process  of  count  in;;,  the  satura- 
tion activity  for  a  ; ivcn  indium  foil  was  obtained. 


o.  «  A.  e  At  (61) 


where  A^  is  the  saturation  act5.\ 

At  Is  t&e  activity  at  I    t, 
^   is  the  decay  constant  for  in' 1 
and   t   is  the  time  after  removal  f*     e  neutron  flux. 

The  saturation  activity  is  proportional  to  the  neutron 
flux  at  an  equilibrium  state.   "he  irradiation  tine  of  the 
indium  foils  for  each  run  was  In  all  cases  at  least  6.5  hours. 
This  length  of  time  corresponds  to  more  than  seven  half-lives 
of  in-.Uum  as  the  half  life  for  indium  is  %k   minutes.  The 
indium  Toils  used  were  not  exactly  all  the  same  size  and 
weight,  rmd   therefore  a  normalizing  correction  was  made.   Li- 
the saturation  activity  for  a  given  foil  by  its  weight 
resulted  in  a  normalized  saturation  activity.    f.nee  a  con- 
stant geometry  was  maintained  in  the  counter,  this  operation 
put  all  of  the  individual  foil  activities  on  the  same  basis. 
This  normalized  saturation  activity  is  the  activity  which  is 
tabulated  in  the  Appendix,  Table  £  through  Table  |0f  for  all 
of  the  foil  positions  in  each  of  the  lattice  confi /"orations 
investigated. 
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Upon  removing  a  foil,  a  three  minute  minimum  decay  tiae 
was  given  to  each  foil  before  beginning  the  count*  This 
period  of  time  eliminated  any  short  lived  decay  components 
in  the  activated  indium,   all  of  the  counts  taken  of  the 
indium  foils  were  three  minutes  in  duration.   The  decay  time 
was  taken  as  the  difference  between  the  time  of  the  foil's 
removal  from  the  neutron  flux  ar     •  mid-time  o£   the  count 
duration. 

The  foil  holder  slot  positions  on  the  front  face  of  the 
subcritical  assembly  are  shown        "e  6.         Uontal 
trection  spacing  was  given  letter  designations   I  rough 
J,   The  vertical  z  direction  spacin  was  numbered  vertically 
1  through  13 •   The  actual  distance  in  in,  from  the  east  face 
was  the  horizontal  y  direction  distance  !csl  latlon.   The 
unit  of  spacing  was  6  in.  in  all  cases  with  the  sole  excep- 
tion of  a  vertical  spacing:  change  to  S»  in.  at  the  position 
where  the  lattice  structure  changed  to  a  5  in.  x  6  in.  area. 

A  glass  wall  Gieger  tube  operating  at  9!>0  volts  was 
used  Tor  all  of  the  foil  activity  determinations.  Before 
and  after  each  series  of  runs,  an  operational  check  on  the 
tube  and  scaler  wa3  made  with  a  sample  of  radioactive  stronti* 
urn.   The  background  activity  remained  essentially  constant  at 
an  average  of     >ounts  per  minute  during  the  series  of  count- 
ing runs.  Dead  time  corrections  were  significant  only  on  a 
few  of  the  higher  counts. 
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The  surveys  taken  o  the  various  lattice  coufiguratlona 
were  as  follows: 
Lattices  I  and  II,  wet  and  dry  cases. 

a.  Vertical  at  x  ■  27  in.;  f  m   13  in.  and  f  m  J 

b.  horizontal,  U   direction,  at  |  ■  18  in.  and 
y  *  30  in.;  r  »30  in. 

c.  Horizontal,  y  direction,  at  x  »  27  in.; 
z  »  JO  in. 

lattice  III,  wet  and  dry  cases. 

a.   Vertical  at  x  »  27  in.;  y  »  JO  in. 
Lattice  without  uranium. 

a.  Vertical  &n&   horizontal,  x  direction,  at 
y  a  3 

b.  Horizontal,  y  direction,  at  x  «  27  in.; 
z   «  30  in. 

Figure  6  is  a  diagram  showing  the  foil  positions. 

\asstone  (5,  .  !  )  shows  t;  .  '  e  fl»  i&tion  in 
the  ve  tical  direction  at  distances  not  too  near  the  top  of 
a  rectangular  suberitlcal  assembly  shape  is 

^(z)    =  ce"yz  (62) 

with 


-ft) 


*2=  (rl2+  (fJ'-O'        i*» 


o 
The  tens  B_  la  designated  aa  the  material  buckli   .   The 

dimensions  of  a  critical  reactor  ean  be  obtained  by  setti 

B_  •  1  c-  the  reometrle  buckling.       the  material  buck* 
■     g 

p 
ling  equal  to  the  geometric  bueklin  ,  B  is  dei'ined  as  the 

bucklin  .   The  base  dimensions  a  and  b  include  the  extrapo- 
lation length.  At  large  distances  from  the  assumed  source 
plane  the  harmonic  corrections  are  assumed  to  be  negligible. 
As  the  neutron  flux  is  proportional  to  the  indium  foil 
saturation  activity,  a  linear  plot  of  the  logarithm  of  fi(z) 
versus  the  vertical  position  z  yielded  the  quantity  ~Y 
which  is  the  slope     riis  plot.   The  base  dimensions  were 
equal  in  the  assembly  used  in  this  investigation.  By  ad 
a:i  extrapolation  1c         .1  A  ^  to  each  side  of  the  l  - 
fectiv©  neutron  "black  wall",  the  square  base     islons  be- 
came  62  in,    qua t ion  63  can  now  be  rewritten  as 

B2  -  2  (A      -  tf2  *  795  x  I)"6  -  *2  (6k) 

where  the  units  of  B  are  cat  . 

mtlon  61;.  was  used  to  evaluate  the  buckling  Tor  the 
wet  and  dry  cases  of  all  lattice  geometries  considered.   The 
slope,  -  X  ,  was  measured  from  the  x  ■  Zl  In*   and  y  «  30  in. 
positions  or  the  vertical  surv>     its  of  each  configuration. 


The  saturation  activities  for  f&Vta  fell  positions  of  the 
configurations  considered  are  tabulated  In  Tables  6  throu-h  10 
of  the  Appendix.   Plots  of  these  values,  Figure  7  thro     i£- 
ure  16,  show  the  vertical  and  horizontal  distribution  of  the 
neutr--     x  Ib  the  subcritleal  aaseribly  for  the  lattice  con- 
urations  Investigated. 

The  theoretical  and  the  experimental  values  of  the  buck- 
ling and  critical  reactor  sl&e  for  the  various  Im 
figurations  considered  are  listed  In  Table  tfe  Figure  17  shows 
the  relationship  between  the  theoretical  and  experimental 
buckling  values  for  the  wet  and  dry  cases  and  the  lattice  | 
Ifig  • 

The  y  »  30  in.  position  of  the  vertical  flux  survey  plots 
was  used  for  the  experimental  determination  I    .oklir     -m 
vertical  flux  distribution  in  the  &   direction  for  the  y  ■  | 
in.  position,     .re  7  and   .'igure  9,  was  plotted  only  to  show 
it  the  flux  variation  at  this  position  was  also  exponential 
in  nature.  Due  to  the  critical  nature  of  t.     ope  of  these 
vertical  flu;-  -lots  with  .respect  to  the  buckling  value,  the 
buckling  from  the  f  m  %     ,  .   ;sitlon  was  not  calculated.   The 
buckli .;  ,   .\uation  6k,   is  a  function  of  the  square  of  the 
slope,  -  X    .  There Tore  a  very  snail  change  in  the  quantity 
- Y   reflects  a  much  larger  change  in  the  value  of  bucklin  . 


, 
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A  deviation  from  a  true  exponential  nature  was  noted  in 
the  vertical  flux  plots  at  the  lower  and  higher  and   positions 
of  all  surveys  in  the  z   direction,      lower  end  deviation 
was  due  to  excess  fast  neutrons,  5  to  12  Kev,  from  the  pin- 
ton!  um-ber;  Ilium  neutron  source.  All  of  the  neutrons 
the  source  in  this  lower  portion  of  the  assembly  had  not  as 
yet  been  moderated.  The  neutron  absorption  reaction  in  indium 
is  for  neutrons  below  1.5'  electron  volts.  The  higher  on^ 
deviation  occurred  approximately  at  the  59  In.  height.    is 
is  th©  point  where  the  lattice  spacing  ofran     or  the  three 
different  ^eometrie  lattices  consider*    I  to  the  construc- 
tion of  the  assembly.  The  slo      the  flux  plots  increased 
near  this  general  position  in  all  of  the  vertical  survey 
plots.   Statistical  error  in  counting  the  small  anount  of 
activity  at  high  z   positions  accounted  for  the  relative  in- 
crease In  the  dispersion  of  the  counts  at  these  positions. 
These  •fleets  can  be  observed  on  the  plotB   shown  on  Figure  7 
through     -e  11. 

The  vertical  flux  i'or   the  wet  conf  it  ur  at  ion  was  greater 
than  the  vertical  flux  for  the  dry  configuration  up  to  a  * 
In.  I'or  f<attice  I,  Figure  ?,  at  the  f  m   Id  in.  position. 
A  similar  situation  existed  for  Lattice  I,   f  ;ro  &a  at  the 
y  «  3    .  position  with  the  cross  over  at  z  m  $Q   i  . 

i  ure  9  shows  that  the  wet-dry  cross  over  point  of  the 
IX  plot  for  Lattice  II  at  y  ■     a.  occurred  at  z  m   30  in. 
Similarly  again,  the  wet-dry  cross  over  point  of  the  vertical 
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flux  distribution  plot  at  y  «  30  In.  Tor  Lattice   ,   -,/ore 
10,  and  Lattice  III,     re  11,  occurred  at  z  «  1+2  la.  and 

k  «  38  In.  respectively.   The  effect  of  the  moderation  of 
the  w&ter  was  greater  than  the  exTeet  of  the  water  absorp- 
tion for  approximately  the  lower  half  of  the  suberltical  as- 
sembly in  all  lattices  Investigated.  The  water  moderation 
affect  increased  the  thermal  neutron  flux  by  slowing  up  the 
fast  neutrons.   In  the  lower  portions  of  the  assembly,  these 
fast  neutrons  are  primarily  from  the  source.  This  moderation 
effect  by  the  water  in  the  process  tubes  was  more  pronounced 
at  the  y  *  30  in.  ositlon  than  at  th©  y  «  18  in*  position 
due  to  the  base  center  position  of  the  source*  Over  the  cr 
fcer,  J   -  30  In.,  the  neutron  flux  consisted  of  more  fast  neu- 
trons at  e  ^Iven  low  horizontal  plane  than  other  positions 
further  off  center*  This  excess  of  fast  neutrons  was  moderat- 
ed In  part  by  th©  water.  At  the  y  «  13  in.  position  a  water 
moderation  effect  was  also  c    '.enced  by  th©  fast  neutrons, 
but  the  greater      ite  mode  ation  due  to  the  greater  dis- 
tance from  the  source  caused  the  number  of  fast  neutrons  to 
be  smaller  than  at  the  eenter  positions*  The  number  ftl    it 
neutrons  due  to  fission  was  also  larger  at  the  center  posi- 
tions than  at  trie  outer  positions  :'oly»  This  is 
due  to  the  sine  dieter    'on  of     thermal  neutron  flux  in  & 
horizontal  plane  in  the  assembly*         neutrons  cau. 
the  fissions  in  th©  natural  uranium  wr-.         -ed  fast  neu- 
trons and  other  fissi          ..    o   position  of  the  mximum 
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thermal  neutr    tax  was  therefore  also  the  position  of  the 
Maximum  fast  neutron  flux  due  to  f issic  . 

The  vertical  flux  curve     ■  the  case  of  the  graphite 
lattice  without  uranium  is  shown  on  Figure  32,  The  dry  con- 
figuration flux  values      attioes  I,  II,  and  III  were  cross 
plotted  to  indicate  the  cross  over  point  on  the  flux  plot 
without  ur     .  The  flux  plot  of  lattice  Hi  crossed  the  no 
uranium  lattice  plot  at  a  slightly  higher  point  than  the  flux 
plot  |    Mice  II    .   'he  flux  of  Lattice  I  had  the  lowest 
cross  over  point  oi'   the  three  lattices*  The  increase  of  neu- 
tron flux  due  to  subcritictl  multiplication  In  Tattle©  I  ap- 
peared to  be  leas  than  the  absorption  effect  of  the  uraniu  . 
This  effect  was  not  nearly  so  pronounced  in  Lat    ■  II  and 

The  neutron  flux  plot  for  the  no  uranium  lattice  ap- 
peared to  be  tmf     :onential  In  character  for  all  but  i 
low  and  high  extreme  positions  In  the  assembly. 

The  slopes  of  the  vertical  flux  plots  for  the  x  *  27  in. 
and  y  »  JO  In,  position  were  determined,  and  it  was  found 
that  the  ma-mitude  of  the;-   slopes  varied  as  follows*   the 
flux  plot  In  lattice  II  had  the  smallest  slope  and  the  flux 
plot  in  Lattice  III  had  the  greatest  slope.   The  slope  of  the 
flux  plot  in  Lattice  I  was  slightly  greater  titan  the  Lattice 
III  slope.  The  greater  the  slop©  of  the  flux  plot,  the  smal- 
ler the  value  o,     buckling  becomes.   This  slope  variation 
was  true  for  both  the  wet  and  dry  configurations  and  was  di- 
rectly reflected  in  the  buckling  values  calculated* 


There  was  little  difference  in  the  effect  of  suberltical 
multipli cation  between  the  dry  cases  of  lattice  II  and  III. 

>m  a  vertical  hi         "  In.  and  upward a,  the  flux  in 
dry  Lattice  II  was  Slightly  I    IV  than  In  dry  Lattice  III. 
The  absorption  differ snet  in  IhtM  two  lattices  was  apparent- 
ly small. 

The  horizontal  surveys  for  eytwaetry  in  the  I  direct!. 
for  x  «  %G   In.  are  plotted  In  Figure  13  *         ire  lf>. 
An  extrapolation  length  was  added  to  each  end   beyond  the  ef- 
fective neutron  "black  wall*.  All  of  the  synsaetry  plots  were 
plotted  as  sine  curves. 

The  horizontal  x  direction  lattice  I  flux,      c  13,  ?or 
the  wet  ease  was  greater  than  the  flux  for  the  dry  case  at 
both  y  »  18  in.  and  y  «  30  in.  positions,  'fhe  horisont&l  x 
direction  Lattice  II  flux,   I  art  13,  showed  that  tnere  was 
very  llttl*         ce  between  the  wet  and  the  dry  cases  for 
either  the  f  m   18  in.  or  y  »  30  In.  positions.  The  y  *  30  In. 
position  flux  for  both  lattices  was  greater  in  magnitude  than 
the  y  »  18  In.  position  flux  due  to  the  sine  distribution  in 
y  direction. 

The  effect  of  the  water  moderation  appeared  to  be  greater 
than  the  water  3b sorption  effect  in  Lattice  I  at  both  the 
y  *  18  in.  and  y  »  30  in.  positions  for  the  x  direction  survey. 
In  the  x  direction  survey  rar   Lattice  II,  these  effects  were 
approximately  equal  for  both  y  positions. 


U6 

Lattice  II  flux  in  the  x   direction  at  z  *  30  in,  had  in 
all  cases  almost  an  ideal  sine  distribution.  The  flux  in 
wet  Lattice  I  had  soot  random  scatter  over-  the  center  section 
of  the  horizontal  survey.     3  variation  perhaps  was  &u®   to 
water  noderatlor;  having  m   g*   ^r  effect  than  the  water  ab- 
sorption over  the  near  canter  section.  This  net  effect  was 
more  pronounced  directly  over  Mm  sour^ 

Plgure  ll|  shows  a  horizontal  x  direction  comparison  at 
z  =*  30  in.  of  the  wet,  dry,  and  without  uranium  eo     •nations* 

,tice  I  values  of  flux  for  both     wet  and  dry  cases  were 
below  the  flux  values  in  the  no  uranium  '.attics.       fan  in 
wet  Lattice  I  was  again  higher  than  the  flux  in  dry  Lattice 
I.   In  Lattice  I,  the  absorption  effect  at  the  %   *  30  ln» 
position  was  obviously  greater  than  the  subcritieal  multi- 
plication effect.   Lattice  II  values  of  flux  for  the  wet  and 
drv  cases  ft*    tst  slightly  below  the  no  uranium  case.  In 

is  lattice  the  subcritieal  multiplication  effect  and  the  ab- 
sorption effects  were  approximately  equal* 

A  lattice  comparison  for  the  horizontal  x  direction  at 
z  ■  30  in.  Is  shown  on  Figure  1  .   The  Lattice  II  flux  was 
greater  than  the  Lattice  I  flux  In  both  the  wet  and  dry  con- 
figurations. 

3y*s?ietry  (sine  distribution)  was  shown  in  the  y  dire 
tion  on  r-if-ure  16  for  Lattices  I  and  II  at  the  x  »  21   in.  and 
z  »  30  in.  position.   The  plot  shows  that  there  was  little 
difference  between  the  wet  and  dry   cases  of  either  Lattice  I 


or  II,      In  Lattice   I  there  was  an  apparent  digression  from 
a  sine  distribution  at  the   center  positions  of   the  wet  con- 
figuration.     Also  the   sine  curve  plotted   for  symmetry  was 
slightly  to  the   left  of  the   points   plotted  on    'igure  16. 
This  was   probably  due   to  the  MMMlt  of   the   Toil  holders 

or  some   spread^  the  tte  blocks   from  a   lattice   I 

loading.      The  flux  plots  of   the  wet  and   dry  configurations 
in  the  y  direction  for   Lattice  II  are  approximately  the  same 
as  the  flux    -.lot  of  the  no  uranium  lattice*      There  Is  apain 
some  dispersion  at  the  center  position   in  thli  arison. 

^he  results  of  the  experimental  and    theoretics!  buckllngs 
determined   for  all  configurations   considered  are  shown  on  F 
ure  17,     The   two  t  tical  mot hods  used  to  determine  the 

buckling  for   the   dry  configurations   showed  reasonable  apre-- 
M&nt,  but  the  wet  case  theory  was   In  disagreement.      The  re- 
sults  from  Method  A  were  hipb   in  comparison  with  the  results 
from  Method  I  for  wet   Lattices    If    II  and   III, 

The  experimental  buckling  values  for  the  wet   cor  a- 

tlons  were  lower  than  the  experimental  buckling  values   for 
the  dry  conripurations  in  all  three  lattice  spacings  Investi- 
gated.     The  theoretical  values   of  buckling  determined  for  the 
dry   conf i {juration st   Methods  A  and  B,   were  a  little  higher 
than  the   experir  |    values   of  buckling  for  the    dry   configure* 

tiona.      The   theoretical  values  of  buckllnp  obtained  ff 
Method   A  for   the  wet  c  ^rations  were  hi^h  in  comparison 

with  the   experimental  buckling   values,      fhe  theoretical   Method 
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B  gave  buckling  values  tor  the  wet  configurations  wh  ch.  were 
In  excellent  agreement  with  the  experimental  buckling  values 
for  the  wet  lattices. 

For  the  natural  uranium  graphite  assembly  s  maximum  did 
occur  for  the  wet  and  dry  cases  of  the  buckling  versus  lattice 
spacing  curve,  T  i^ure  17.   This  maximum  occurred  in  both  the 
experimental  and  theoretical  analyses*   The  C/tf  aton  ratios 
were  73  for  Lattice  I.  ll|£  for  Lattice  II,  and  292  for  Lattice 
III.  Davenport  (8,  p.  3155  concluded  that  the  maximum  buck- 
line  for  a  given  natural  uranium  fuel  assembly  should  occur 
at  a  value  of  c/U  between  £0  and  100*  The  maximum  buckling 
in  this  investigation  of  three  lattice  spacings  occurred  at 
about  a  C/U  value  of  Ifcft  for  both  the  wet  and  dry   c     ura~ 
tlons. 
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VI.   DI3CU3SI0N  OF  1     73 

The  theoretical  methods  employed  had  the  following  prob- 
able errors.   Accurate  detailed  neutron  scattering  and  capture 
cross  section  data  are  not  available  in  the  literature  for  all 
materials.   Corrections  for  varying  densities  and  purities  of 
the  materials  were  not  computed.   A  calculated  disadvantage 
factor  for  uranium  was  used  instead  of  a  measured  value.  Val- 
ues from  the  literature  were  used  for  the  ,'ermi  age  and  the 
various  lattice  diffusion  lengths  where  experimental  values 
again  would  have  given  more  accurate  results  for  the  calcu- 
lated buckling.  There  were  cracks  between  the  graphite  blocks 
with  a  process  tube  uranium  load  in  the  assembly.   These  air 
spaces  and  the  air  effect  around  the  process  tube  itself  were 
assumed  to  be  negligible.  The  effect  of  temperature  and  hu- 
midity on  the  results  was  also  assumed  to  be  negligible.  Fi- 
nally, one  f-roup  neutron  theory  was  used  throughout  the  analy- 
sis* 

The  following  comments  can  be  made  cone «  rnlng  the  experi- 
mental methods  employed.   The  harmonic  corrections  to  the  flux 
distributions  obtained  were  assumed  to  be  negligible.   The 
positioning  of  the  foils  and  the  normalizing  of  the  foil  ac- 
tivity were  subject  to  human  error,   r inally.  and  most  impor- 
tant, the  slope  of  the  In  ft[z)   versus  2  direction  curve  must 
be  determined  with  great  accuracy.   This  slope  value  is  very 
critical  in  determining  the  actual  buckling.   Repeat  runs  for 


consistency  at  a  given  position  should  be  taken  from  which 
an  average  value  of  activities  at  a  position  could  be  used 

for  plottinr  ard  the  buckling  calculation*  Varying  the  Toil 

4 

weight  at  a  giv*n  position  during  the  runs  would  tend  to  can- 
cel out  any  normalizing  error* 

The  feasibility  of  an  actual  reactor  with  the  process 
tubinr,  natural  uranium,  slug  geometry,  and  water  coolant  uaed 
in  this  investigation  was  determined.   Assuming  1250  channels, 
a  water  coolant,  and  0*^05  in,   flow  area  per  channel  the  mass 
flow  is  1*68  x  10  lb  /hr.  Assuming  a  temperature  difference 
of  70*P  and  using  C  ■  1.06  for  water  at  1^00 *F,  the  heat  rate 
transfer  is  366  megawatts.   This  heat  rate  transfer  is  large 
enough  for  a  potential  power  reactor*   The  above  assumptions 
are  conservative  In  comparison  with  the  experimental  critical 
reactor  volume  determined  for  the  wet  Q*5   in*  lattice  config- 
uration. 

The  process  tubing  did  spread  the  graphite  blocks  a  small 
aitount  due  to  the  fact  that  some  of  the  holes  were  smaller 
than  the  normal  average  designed  size.   A  power  drill  of  an 
appropriate  size  could  enlarge  these  holes  so  that  the  process 
tubing  would  slide  in  and  out  of  the  assembly  with  ease  and 
no  spread  would  occur.   It  was  found  that  a  lattice  change  was 
more  easily  accomplished  using  the  process  tube  uranium  con- 
tainment than  with  Just  the  loose  slugs.   One  end   of  the  proc- 
ess tube  was  beveled  to  facilitate  loading  of  the  wire  wrapped 
slugs.   It  is  suggested  that  a  number  seven  rubber  stopper  be 


51 

well  seated  In  the  beveled  end  and  a  larger  stopper  be  need 
at  the  other  end.   Coolant  leakage  would  in  thla  way  be  elimi- 
nated. 

The  E-3  posiMon  foil  holder  hole  was  a  through  hole  in 
the  y  direction,  This  through  hole  would  be  of  wore  use  if 
it  were  moved  up  to  the  E-8  or  E-9  positions*   The  E-3  posi- 
tion was  too  close  to  the  source  for  a  good  symmetry  flux 
survey. 

The  8«5>  in.  lattice  is  the  optimum  lattice  spacing  which 
is  possible  with  the  subcritical  assembly  used.   Suggested 
Ideas  for  future  study  are  the  investigation  of  other  coolants 
and  a  more  detailed  analysis  of  the  S.£  in.  lattice  for  both 
the  wet  and  dry  configurations.   Also,  the  effect  of  the  re- 
flection by  the  water  contained  below  the  source  may  be  in- 
vestigated to  rive  a  maximum  source  flux  with  the  available 
sources  in  place.   With  the  use  of  an  additional  inner  process 
tube  the  effect  of  a  variable  water  annulus  could  also  be  in- 
vestigated. 


octe* 
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VII.         IONS 

The  moderation  effect  of  the  water  was  greater  than  the 
absorption  effect  of  the  water  for  approximately  the  lower 
half  of  the  autcritlcal  assembly  for  the  three  lattice  spac- 
ings  investigated.   The  thermal  neutron  flux  in  the  lower 
portion  o£   the  assc.i&ly  was  greater  in  the  wet  cc     rations 
than  in  the  dry  configurations,  and  the  reverse  situation 
existed  in  the  upper  portion  of  the  mssc-mbl    The  6  in.  wet 
lattice  had  the  lowest  overall  flux  because  the  absorption 
effect  of  the  greater  amount  of  water  and  uranium  was  larger 
than  the  subcritical  multiplication  effect  in  this  lattice. 
The  12  in.  lattice  was  *over  moderated"1  with  an  apparent  snail 
•uberitical  multiplication  effect.   The  optimum  lattice  in- 
vestigated was  the  $*$   in.  lattice.   In  this  lattice  the  water 
had  very  little  effect  on  the  flux  distribution  due  to  the 
balance  of  the  water  moderation  and  absorption  effects. 

The  aysnwtry  surveys  showed  that  in  the  x  and  y  direc- 
tions that  the  flux  did  follow  a  sine  distribution  for  both 
the  wet  and  dry  configurations. 

The  two  theoretical  methods  used  to  determine  the  buck- 
ling for  the  dry  configurations  gave  results  with  reasonable 
agreement.  These  theoretical  buckling a  for  the  dry  cases  were 
high  in  comparison  with  the  experimental  bucklings  determined 
for  the  dry  cases  except  for  the  12  in.  lattice.   The  two 
theoretical  methods  used  to  determine  the  buckling  for  the 
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wet  configurations  gave  results  which  were  in  disagreement. 
The  method  of  Murray  gave  values  or  the  buckling  for  the  wet 
configuration  of  the  6  in,  and  8.5  in. lattices  which  were 
very  high  in  comparison  with  the  experimental  buckling  values* 
The  method  of  Rumaey  and  Volkof t   gave  values  of  the  buck- 
ling for  the  wet  configuration  of  the  6  in.  and  8.5  in.  lat- 
tices which  were  in  excellent  agreement  with  the  experimental 
buckling  values  determined.   The  experimental  buckling  value 
determined  for  the  12  in.  wet  lattice  was  between  the  two 
theoretical  buckling  values  determined  for  this  lattice.  The 
experimental  buckling  values  for  the  dry  oases  were  higher 
than  the  values  for  the  wet  cases  in  all  lattice  spacing* 
investigated.   The  only  experimental  positive  value  of  buck- 
ling for  the  wet  configurations  was  the  buckling  value  for 
the  3.5  in.  lattice. 

A  maximum  in  the  buckling  versus  lattice  spacing  curve 
did  occur  for  both  the  wet  and  dry  configurations.   This 
amxisRim  in  buckling  was  found  in  both  the  theoretical  and  the 
experimental  values  determined.   These  buckling  curves  for 
the  wet  and  dry  configurations  did  not  cross;  hence  a  "fail- 
safe" behavior  upon  the  loss  of  water  coolant  could  not  be 
expected  from  these  results. 
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Table  1 
Dimensions  for  a  unit  cell 


lattice  I 

(6  in.),  II  (8.5  in.),  and  III  (12  in. 

) 

uranium  rod  radius 

1.270 

M 

fc.i » 

thickness  of  aluminum  I 

can 

. 

M 

*w> 

effective  thickness  of 

water  annulus      .273 

V 

effective  thickness  of 

process  tube 

.102 

cm 

'air  ' 

r  effective  thickness  of 

air 

annulus 

o.k5S 

en 

rl  ' 

equivalent  inner  radius  ox 

graphite 

2.. 

cm 

iv&lent 

cross  section  areas  for 

one 

cell 

5.o6 

Uranium 

cm2 

Slug  can 

.3U 

2 

M 

Water 

ooolant 

2.61 

2 

cm 

Process  tube 

.06 

cm 

Air 

5.65 

cm 

Graphite,  Lattice  I 

m 

217. 

cm 

ophite,  Lattice  II 

Ufc 

Mr 

Graphite,  Lattice  III 

66C. 

cm 

Volumj  •■■■,;■ 

Vu  • 

1*0 

cmr 

v«l  • 

slug  can  and  cap 

23.0 

cm** 

Vw  • 

water 

55.7 

CUT 

VP- 

process  tube 

22.6 

CRT 

v 

gl  ' 

aphite,  r^ttlce  I 

25 

V  ■ 

graphite,  Lattice  II 

9250 

mr 

V«3  - 

grap;    ,  Lattice  I 

1G,500 

cur 

'• 
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Table  1  (Continued) 


Lattice  Xf  r2,  equivalent  outer  radius  of  graphite  8.60  cm 

Lattice  II,  r?,  equivalent  outer  radius  of  graphite  11.9L}.  era 

Lattice  III,  r2,  equivalent  outer  radius  of 

graphite  16.76  era 


Symbol 


Table  2 
Lattice  cell  material  constants 


Value 


Source 


I: 

I? 

It 


tf 


u 


r 


0.01323  e» 


-1 


0.01556  cm 


-1 


i.32h  cm' 


-1 


0.399  cm 


-1 


0.0967  cm 


-1 


-1 


0.00036   cm 
0.375  em"1 

.   1057  cm"1 
0.35U  cra"1 
0.017  era"1 
0.675  cm"1 


0.01992 


-1 


(5)    and  liquation  25 
(5),    (ID    and   Equation 
(1)    and   Equation  2$ 
(1)    and  Equation  25 

uation  \& 
(5) 

(1)  and  Equation  2$ 
(1)  and  I  qua t Ion  l^. 
(1)  and  Equ&tion  ij.9 

(5) 

Equation  26 

(5)  and  F qua t ion  27 
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Table  2   (    ontinued) 


Symbol 

Value 

Souroe 

*w 

•  3U72  em"1 

(5)   and  Equation  27 

0.U218  cm"1 

Equation  I4.6 

*r 

).1059  oaf1 
0.0555 

Equation  1)7 

COB    £ 

-uation  50 

«A 

20 

(3) 

I 

%158 

(5) 

Fermi  age. 

r     life  era2 

(1) 

4 

g 

2520  cm2 

(1) 

fi 

1*56  g/csr 

Measured  value 

Cu 

19.0  g/cm3 

Measured  value 

3/fc 

.0828  em2/gm 

alculated  value,   Equation  \\Z 
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-e 


'"*  | 


■) 


<9B 
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CO 
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Table  3 
Theoretical  nuclear  lattice  conatanta 


Method 

A 

I 

II 

0.926 

.370 

0.823 

0.909 

1.023 

1.062 

0.901 

0.8^7 

0.823 

0.909 

0.996 

i.)35 

Method 

. 

Lattice 

DRY 

t  .926  0.870                              .767 

P  0.823  0.909                              .955 

kp.*  1.023  1#062 

f  .  >oi  0,81^7                           .751 

p  0.823  0.909                              .955 


Lattice  I  II  III 

_ 
f 

P 


-.920 

0.863 

0.759 

0.836 

0.915 

.958 

1.03U 

1.061 

.976 

0.882 

.  318 

.710 

.  .32 

0.911 

•fSk 

0.985 

1.008 

. 

.luated  us:     6  m   1.  )27   and  Y\  *  l. 


7 


or, 
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Table  k 
critical  reactor  size 


Theoretical 

Method  A 

rjattic© 

I 

II 

ill 

DRY  Value 

2cm2 
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1*1  x  mT* 

m  x  icf6 

-17  s  IO"6 
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611  x  106 

ft   106 
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19.2 

■ 
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r2     2 
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•       e       **»**• 

Theoretical 

26. k 
He t hod  B 
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I 

II 
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63  x  io"6 
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»  eft"*™»ft 
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MM 
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ft. 
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19.6 

mm 
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ij  ess 
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kS9 
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2     -2 
B^  em 

-zy  x  10~6 
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-  H  x  io"6 
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M 

53  ko  x  106 

« 
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Table  l^  (Continued) 
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Lattice 
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Figure  3*   Suboritlc&l  assembly 
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Figure  7.   Vertical  flux  survey,  Lattice  I, 
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